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Summary 

A knowledge of building wake characteristics is useful in a variety of applications 
including dispersion of pollutants downwind of conventional or nuclear power plants, 
airport runway interference effect, take-off/landing limitations at heliports, pedestrian 
wind comfort, and wind loads on structures. Recent investigations have advanced both 
understanding of physical flow processes occurring in the near and far wake regions 
and theoretical predictive capability for flow and diffusion in the far wake region. These 
studies have shown, for example, that separation cavities immediately downwind are not 
enclosed by free streamlines as in two-dimensional separation/reattachment, that organized 
vortices can play an important role in wake development, and that theoretical solutions 
for some wake regions are possible. 


Notation 


A reference area for concentra¬ 

tion measurements 
a virtual origin of wake 

C Fx force coefficient of building 

in x direction 

C w wake concentration of pol¬ 

lutant 

H building height 

h vertical height of center of 

vortex 

K von Karman constant 

k 2 constant 

n mean velocity power law 

exponent 

Q pollutant source strength 

T mean fluid temperature 


T H mean temperature at height 

H in approach flow 
T Zi mean temperature at the 

height Zj 

U mean velocity in x direction 

U H approach mean velocity at 

height H 

U b approach mean velocity at 

height 6 

u mean velocity defect in 

wake (U(z) approach -U(z) 
wake) 

x,y,z space coordinates, x down¬ 

wind, y lateral, z vertical 
y v lateral distance to center of 

vortex 
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Greek symbols 

a approach wind direction 
r circulation of vortex at jc = 0 

5 boundary-layer thickness 


6 mean temperature excess in 
wake ( T(z) wake - T(z) ap¬ 
proach) 

7 , X constant of order one 


Introduction 

During the past decade, much has been learned about the characteristics 
of flow about buildings and the structure of the wake downwind of buildings. 
Knowledge of building wake characteristics is useful in a variety of applica¬ 
tions including dispersion of pollutants downwind of buildings, such as: con¬ 
ventional or nuclear power plants; airport runway interference effects; take¬ 
off/landing limitations at heliports; pedestrian wind comfort; and wind loads 
on structures. Recent investigations have advanced both the understanding 
of physical flow processes occurring in the near and far wake regions and the 
theoretical predictive capability for flow and diffusion in the far-wake region. 

There is a relatively high level of misunderstanding of how wind flows 
around buildings, probably caused by conceptual extensions of two-dimen¬ 
sional flows. For example, Fig. 1 shows that the separation zones for a two- 
dimensional obstacle are bounded by streamlines so that the separated cavi¬ 
ties are closed. Thus, pollutants released in the separated cavities can escape 
only by turbulent diffusion across the cavity boundary. For a three-dimen¬ 
sional object, the separation lines of Fig. 1 are not correct as will be dis¬ 
cussed in more detail below, and pollutants may enter or escape the separa¬ 
tion cavity by convection. Some prediction formulae incorrectly assume that 
pollutants can enter the three-dimensional separation cavity only by turbu¬ 
lent diffusion across this cavity boundary. 


CORRECT: 2-0 BODIES 
INCORRECT: 3-0 BOOIES 



Fig. 1. Separation cavities for two-dimensional flow. 


Near-field flow 

Very little data defining the flow field about buildings are available from 
full-scale tests: Colmer [1] and Frost and Shahabi [2] reported limited 
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measurements of mean velocity and turbulence intensity downwind of 
buildings from which decay of velocity defect with distance downwind could 
be determined. Other full-scale measurements related to pedestrian comfort 
have reported velocities at isolated points near buildings, but not in a way 
which would contribute significantly to an understanding of the flow char¬ 
acteristics about buildings. 



Fig. 2. Surface shear-stress patterns [8,9], 
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Much of the current knowledge relating to three-dimensional wind flows 
around buildings has been obtained by means of model tests in boundary- 
layer wind tunnels. Important publications based on such wind-tunnel tests 
are refs. 3—10, and others too numerous to include in this brief review. An 
excellent review including more details related to diffusion is presented by 
Hosker [11]. These studies showed that three-dimensional obstacles cause 
fundamental differences in flow characteristics from those originated by 
two-dimensional obstacles. Oil flow studies [8, 9] (see Fig. 2) permitted an 
identification of mean surface-shear-stress patterns on both the ground plane 
and on the building surface. In combination with smoke flow studies and 
kinematic analysis, these data provided the basis for the mean flow patterns 
which are shown in Figs. 3—5 based on refs. 8 and 9. 



Fig. 3. Mean streamline patterns about a building [8,9 ]. 



Fig. 4. Centerline streamline patterns for flow reattaching to top [9], 



Fig. 5. Centerline streamline patterns for flow not reattaching to top [9], 
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Figure 3 shows that the mean streamlines for flow about a single building 
form a complex flow geometry. The flow approaching the obstacle separates 
from the surface at a distance upwind of the building. This separation loca¬ 
tion is dependent, to the first order, on building height-to-width ratio, 
building height-to-boundary-layer-height ratio and upstream surface rough¬ 
ness. The vorticity in this separated flow, in combination with the pressure 
distribution on the front of the building, results in downward flow on the 
front of the building. This subsequently causes the separated flow to roll up 
into a vortex. The vortex is wrapped around the building by convection into 
a horseshoe shape. The primary vortex induces additional vortices of smaller 
size and strength (Figs. 3 and 4) which are eventually incorporated into the 
primary vortex around the side of the building and lose their individual iden¬ 
tity. The horseshoe vortex can be identified in the flow at some distance 
downwind. 

The wind that impinges on the front of the building forms a stagnation 
region between 2/3 and 3/4 of the building height depending on building 
height-to-width ratio. From this region, flow moves outward toward all front 
edges of the building. The flow separates at the front edge of the top and 
sides and may or may not reattach to the top or sides before reaching the 
back edges. Reattachment depends on building length-to-width ratio, height- 
to-length ratio and upstream roughness (which determines the turbulence 
intensity in the approaching wind — a significant factor in distance to re¬ 
attachment). Figures 3 and 4 show flow patterns for a reattached flow; Fig. 
5 shows flow patterns for an unreattached flow. 

A separation cavity covers the rear face of the building and also the top 
and sides when reattachment on the sides does not occur. Examination of 
Figs. 3—5 reveals that the mean streamline reattaching to the ground surface 
downwind of the building originates in the flow upwind of the building and 
not at a separation point on the building. The same is true of the streamline 
which forms a stagnation point on the rear face. Thus, the mean flow field 
convects mass directly into the separation cavity, and, by continuity of mass, 
it also convects flow out of the cavity region. The direct convection out is 
more difficult to visualize but involves off-centerline flow departing from the 
two standing vortices with vertical axes forming an arch at the back corner 
and probably entrainment into the horseshoe vortex. For this reason, trans¬ 
port of pollutants into or out of the separation cavity is primarily accom¬ 
plished by convection and not by turbulent dispersion. If the cavity length is 
defined by the centerline reattachment point downwind of the body, then 
the cavity length can vary from 2 to 6 or more buildings heights [11]. Be¬ 
cause of the direct convection into and out of the cavity, a precise definition 
of cavity boundary may not be meaningful. 

The organized picture of the cavity described above is made cloudy by 
the existence of a high level of turbulence in this region. Smoke flow and 
bubble track experiments suggest that the separation cavity is completely 
“washed out” by turbulent gusts at times. Thus, the use at this time of the 
mean streamline model described above for more than a qualitative under- 



26 


standing of the flow field must be approached with caution. Some quantifi¬ 
cation of the mean flow field has begun with velocity-measuring devices 
capable of correctly evaluating the effects of high turbulence levels and 
reverse flows (see ref. 12). 

The flow field described in Fig. 3 has significant implications for pedestri¬ 
ans and for the dispersion of pollutants about the building. For example, the 
wind flows downward on the front face of the building from the front stag¬ 
nation region and around the comer of the building at the base as it wraps 
up into the horseshoe vortex. This effectively brings velocity magnitudes 
associated with approach wind velocities at about 3/4 the height of the 
building to the surface at the front comers of the building. Immediately 
around the comer at ground level in the separated flow zone, the velocities 
are quite low. Most city pedestrians who have walked around the corner of a 
tall building into a blast of wind are familiar with this phenomenon, if not 
its cause. The addition of nearby buildings can, in some cases, increase this 
velocity. The high turbulence levels associated with the separated zones can 
cause a decrease in control for helicopters which might try to land on the 
roof or ground. Pollutants from roof exhausts can become entrained in flows 
which recirculate to the back and sides of the building where ventilation in¬ 
takes are frequently located; see, for example, ref. 13. A significant amount 
of indoor air pollution and cooling-tower recirculation can be attributed to 
the circulations within separated flow zones about buildings. 

Far-field flow 

The reattached flow downstream from the building requires some distance 
before the wind recovers the characteristics of the approach boundary layer 
and all disturbances caused by the building disappear. This wake region is 
characterized when compared to the approaching wind flow by a lower mean 
velocity, higher turbulence intensity, smaller scales of turbulence and, in 
some cases, organized vorticity with its axis parallel to the wind direction. 
Figure 6 defines several variables used below. 



Fig. 6. Definition of variables. 
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Figure 7 shows the distribution of velocity defect (defined as the differ¬ 
ence between upstream approach velocity at height z and the local mean 
velocity at height z) on the building centerline downwind of a building from 
ref. 5 with H/S = 0.11, H/W = 0.4, W/D = 3.25 and a = 0°. The maximum 
velocity defect gradually decreases but is still detectable at x/H = 18 for this 
combination of building and approach flow characteristics (power law expo¬ 
nent = 0.24). A slight increase in velocity is noted just above the separated 
shear layer forming the top edge of the cavity for x/H = 0.5—2. The wake 
height increases monotonically with distance downwind to a height of 3—4 H. 


Y= 0 



0 0 0 0 0 0000000 


0 0.2 0.4 u a 

Fig. 7. Vertical profiles of mean velocity defect on wake centerline [5], 
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Y/(W/2) 

Fig. 8. Horizontal profiles of mean velocity defect in wake [5], 
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Figure 8 (from ref. 5) shows the distribution of mean velocity in the hori¬ 
zontal direction perpendicular to the approach velocity at a height of zjH = 
0.8. The same characteristics can be observed which were noted for Fig. 7: 
a gradual decrease in magnitude of velocity defect; an accelerated flow near 
the separated shear layer; and a monotonically increasing wake width to 
about 4W/2 (four half-widths of the building from the centerline). A second¬ 
ary maximum in the velocity defect profile can be noted when xjH is from 2 
to 12 near y/{WI2) = 2—2.5. This is the result of the horseshoe vortex which 
convects lower momentum fluid upward in the boundary layer at the region 
of the secondary maxima increasing the velocity defect. The vortex also con¬ 
vects higher momentum fluid downward in the region between the maxima 
near y/(W/2) = 1—1.5 resulting in a smaller velocity defect at that location. 

The wake characteristics of Figs. 7 and 8 are not universal for all building 
aspect ratios or approach roughness conditions. A rougher boundary upwind 
of the building results in a higher level of turbulence in the approach wind 
causing the wake to disappear more quickly and evidence of the horseshoe 
vortex in the wake to diminish. A smoother approach boundary causes the 
wake to be detectable to a greater downwind distance and the horseshoe 
vortex to remain in evidence longer. 

The influence of building aspect ratio on magnitude and persistence of the 
wake is presented in Fig. 9, modified from ref. 5, in which the largest velocity 
defect at any downstream distance is plotted against distance downwind. 
The model data were acquired in a 0.24 power-law boundary layer. Field data 
from Colmer [1] and Frost and Shahabi [2] are plotted with the wind- 
tunnel data. For practical purposes, the extent of the wake in a suburban 
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Fig. 9. Maximum mean velocity defects in the wakes of several buildings [5 ]. 
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setting is 3—6 building heights for buildings with width 1/3 or 1/4 their 
height and is 10—30 building heights for buildings whose width is 1—4 times 
their height. The Colmer data may not have been obtained on the wake center- 
line. This would result in an underestimation of the wake defect. 

One characteristic of building wakes is usually a substantial increase in 
turbulence in the wake, accompanied by a shift in turbulent energy to 
smaller scales. Figure 10 shows the decay of turbulence downstream from 
the model buildings whose mean velocity defects were plotted in Fig. 9. The 
effective length of the turbulence wakes was about the same as for the mean 
velocity wake. The turbulence increase in the wake downstream is significant¬ 
ly reduced when the vortex wake is stronger; for example, in the wake of a 
hemisphere [6] or in stably stratified flow [10]. In this case high-momen¬ 
tum, low-turbulence flow from above the hemisphere is brought down close 
to the ground on the wake centerline by the vortex pair and results in a 
mean velocity excess and a turbulence defect. 



X/H 

Fig. 10. Maximum turbulence intensity excess in the wakes of several buildings. 

A rectangular building can produce vorticity in the wake with axis parallel 
to the approach flow direction by mechanisms other than those associated 
with the horseshoe vortex. With a building rotated so that one side is approx¬ 
imately 45° to the approach flow, a pair of rooftop vortices forms over the 
upwind comer whose axes can become aligned with the flow in the wake. 
When one of these vortices is much stronger than the other, its presence in 
the wake can dramatically extend the length of the wake. Figure 11, from 
ref. 14, shows the mean velocity defect downwind of a building whose side 
ratio is 2.6. A mean velocity excess remained at a significant magnitude 80 
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Fig. 11. Horizontal profiles of velocity defect in the wake of a building at angle to the 
flow [14]. 

building heights downwind where wind-tunnel test-section length limitations 
prevented further measurements. 

Theoretical predictions 

Despite the considerable complexity of the flow around and downstream 
of a building, a reasonably simple but effective analysis is possible to provide 
quantitative predictions of wake behavior. In 1969, Hunt and Smith [15] 
published a momentum-wake theory for the mean velocity defect in the 
wake of a three-dimensional obstacle immersed within a turbulent boundary- 
layer flow. Momentum wake theory is based on the assumptions that only 
small perturbations to the turbulent boundary layer are caused by the ob¬ 
stacle, that the ratio of the model height to the boundary-layer height is very 
small, that the velocity profile power-law exponent is small, and that the 
mixing may be represented by a constant eddy viscosity. This theory is valid 
only for the far-wake region because the separated region behind the 
building cannot be considered to be a small perturbation on the approach 
flow. The mean velocity defect predicted from the three-dimensional 
momentum-wake theory can be expressed as (see Fig. 6 for coordinates and 
wake schematic) 

u(x ’. y ' z) = k l F 1 (z",y")/[(x - 

Uh 


( 1 ) 
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where u(x,y,z) = velocity defect, U H = approach velocity at building height 
H, a = virtual origin of the wake, n = velocity profile power-law exponent, k : 
is given by 

fe, = 0.21 C Fx (W/H)/{4:X V2 y[2K 2 n/(2+n)] li3+n)K2+n)] } 

where W = width of building, X and y ~ constants of order one, K = von Kar- 
man constant (0.41), C Fx = force coefficient of building, and Fx(z",y") is 
given by 

„ , „ if, 7? [1/(2+n)I exp{- [r? + y" 2 /(1.5 + r?)] } 

F ' ( * • y » * -hTTip- 


where 


z" = (z/H)/[(x-a)/H] llK2+n)] 

n = z"< 2+n >/[2(2 + n)K 2 ny] 

" = (Wg) / 2 + n 

~ V _ [(x - a)/H] t 1/ < 2+n)1 ' 2XyK 2 n ' 


Thus, for the wake of a three-dimensional building and small n, the theory 
predicts that the mean velocity defect decays as x" 1 ' 5 . A similar prediction 
for the turbulent intensity excess decay rate suggests proportionality to x' 2 
[16]. 

Measurements of Peterka and Cermak [5] did not confirm the existence 
of a Gaussian distribution of velocity defect in the crosswind direction in¬ 
corporated in the momentum-wake theory. They concluded that this effect 
was the result of the horseshoe vortex which forms in the stagnation region 
at the front of the building and wraps around the building to produce a 
mean vorticity with the axis in the longitudinal direction. Noting this result, 
Hunt [17] subsequently developed a vortex-wake theory to account for 
these vortices. The longitudinal velocity excess in the inviscid zone was 
determined by Hunt [17] and given by 


u 


nThy’x 

TT 


' _1_ 

-(y >2 + h 2 + z 2 ) 2 - Az 2 h 



( 2 ) 


where y' = y v - y. y v and h are the distances to the center of the imposed 
vortex in the lateral and vertical directions, respectively, and V is the circula¬ 
tion at x = 0. 

The physics of vortex- and momentum-wake mechanisms were applied to 
any passive scalar such as temperature by Kothari et al. [10], They devel¬ 
oped the theory to predict the time mean temperature field behind a three- 
dimensional bluff body immersed deeply in a turbulent stably stratified 
boundary layer. Some of the results for the velocity field from Hunt’s 
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vortex-wake theory [17] were used to derive the model for the temperature 
field. 

Transport of the scalar quantity temperature is governed by the energy 
equation. With the assumptions that only small perturbations in the thermal 
boundary layer are caused by an obstacle, that the ratio of the model height 
to the boundary-layer height (velocity or temperature) is very small, that the 
velocity and temperature power-law exponents are small, that the dissipation 
and molecular diffusivity are small, and that the mixing may be represented 
by a constant eddy diffusivity, Kothari et al. [10] obtained the temperature 
perturbation induced by the obstacle on the approach thermal boundary 
layer. The three-dimensional temperature perturbation solution for a flow 
field subjected to only momentum-wake effects is given by 

e(x,y,z) = T h F 2 (z", y") 

DELTA T 2 DELTA T / x - a >/<»*» 


where DELTA T = T s - T z in the undisturbed boundary layer, T 6 = mean 
temperature at the top of the turbulent boundary layer, T z = mean tempera¬ 
ture at a reference height z t , T H = mean temperature at height H in the 
approach flow, k 2 = constant determined from experiment, a - virtual origin 
of the wake, n = velocity profile power-law exponent, F 2 (z" ,y") is given by 


F 2 (z",y") = 


,[l/(2+rc)] 


(V + 1.5) 1 


(V + 1-5) 


where 



z "(2+n) 

[2(2 + n)K 2 n0.85y] 


„ (y/H) r (2 + n) -| l,/(2+n>1 

y ( x -q \ [1/(2+n)] L2\0.857lC 2 n J 
l B ) 

K = von Karman constant (0.41), c 8 = constant (not needed at y - 0). 

The excess temperature induced by the trailing horseshoe vortices was also 
determined by Kothari et al. [10] and is given by 
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To consider the effect of two horseshoe vortices along the centerline of a 
building the perturbation effects of the vortices must be doubled in the cal¬ 
culation. The total velocity defect and temperature excess are given by the 
algebraic sum of eqn. ( 1 ) and twice eqn. ( 2 ) and by the algebraic sum of eqn. 
(3) and twice eqn. (4), respectively. 

Comparison of measurements with theories 


The theories of Hunt [15,17] and Kothari et al. [10] for the mean veloc¬ 
ity and temperature for momentum and vortex wakes behind the building 
are evaluated in this section. The experimental data of Peterka and Cermak 
[5] and Kothari et al. [10] have been utilized for comparison with the 
theories. 

The momentum theory requires knowledge of the force coefficient, C Fx , 
on the building which was derived from data of Akins et al. [18], The con¬ 
stants X and 7 were assumed equal to one. The virtual origins for the build¬ 
ings were derived from the technique suggested by Hunt [16]. 

From the momentum-wake theory, eqn. (1), it can be shown that the 
maximum velocity defect occurs on the wake centerline at 


z 


H 



(5) 


The maximum velocity defect from eqn. (1) at this height for small n and 
large x is 


u W/xY i/2 f X(H\ 2/ *- 

u H 01 h\h) ~ W/ . 


( 6 ) 


The constant of proportionality can be obtained by evaluation of eqn. (1). 
The authors have found that an additional factor to account for initial verti¬ 
cal mixing in the near wake improves the correlation of data obtained behind 
rectangular buildings 



(7) 


Data from Fig. 9 are plotted in Fig. 12 using variables defined by eqn. (7). 
Also plotted in the same variables is the range of maximum velocity defects 
predicted by eqn, (1) using eqn. (5) and C Fx = 1.2, 7 = X = 1, n = 0.2 and a = 
-2H. The data collapse well, but fall under the momentum wake predictions. 
The difference between prediction and data may be due to vortex-wake 
effects which could be predicted using eqn. ( 2 ) or due to initial vertical 
mixing in the near wake not included explicitly in the theory. Application of 
eqn. ( 2 ) did not give a consistent improvement for all cases. Until further 
study is completed, a reasonable approximation of peak momentum defect 
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for rectangular buildings is 


u 

Uit 



1.3 


in which a constant full-scale value of 5 = 1000 ft. (~300 m) can be used 
without significant impact on the answer computed. 



Fig. 12. Comparison of maximum momentum defect with theory. 


The data of Kothari et al. [10] in a stably stratified flow showed indica¬ 
tions that the horseshoe vortex was an important factor in the wake struc¬ 
ture. Figures 13 and 14, from ref. 10, show the comparison of analytical and 
experimental vertical profiles of mean velocity and temperature on the wake 
centerline. The calculations of momentum- and vortex-wake effects were 
made using the theories discussed above. The vortex location was estimated 
from data in ref. 6, while vortex strength was estimated by a procedure 
developed in ref. 10. The effect of the momentum wake is to decrease veloc¬ 
ity and decrease temperature while the effect of the vortex wake is to in¬ 
crease velocity and temperature by pulling high-momentum, high-tempera¬ 
ture fluid down into the wake. As shown in Figs. 13 and 14, the sum of the 
momentum- and vortex-wake effects matches the data well. One constant in 
the temperature theory, ki, was not available through the theory and had to 




35 


be set by comparison with data. The value of k 2 was selected to force a 
match between theory and experiment at x/H =10 and z/H = 1. Additional 
research is needed to determine a method for predicting k 2 without data. 
No free constants were needed for the velocity comparisons. The comparison 
between theory and experiment is excellent except in the near wake where 
the perturbations are large and the theory should not be valid. Note that the 
vortex is sufficiently strong in this weakly stably stratified flow that a veloc¬ 
ity excess and temperature increase are produced in the wake. The momen¬ 
tum-wake theory alone cannot adequately model such a wake. 



Uh 


Fig. 13. Comparison of vertical profiles of velocity defect on wake centerline in stably 
stratified flow with theory [10]. 



Fig. 14. Comparison of vertical profiles of temperature excess on wake centerline in 
stably stratified flow with theory [10], 
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Concentration prediction 

A theory for dispersion of pollutants in the wake has been developed [19] 
in a manner parallel to that for temperature discussed above. Application of 
this theory to dispersion in the wake of nuclear power plants has been made 
[20,21], Figure 15, from ref. 20, shows a comparison of largest mean con¬ 
centration in the wake of a power plant and the prediction. One free coeffi¬ 
cient appears in the theory. Comparisons between theory and experiment to 
date are encouraging and are the subject of continuing research. 



Fig. 15. Comparison of ground-level concentration on wake centerline with dispersion 
theory [20]. 

Conclusion 

On the basis of research conducted during the last decade investigating 
wind flow around and downwind of buildings, several conclusions are evi¬ 
dent. 

(1) Wind flow about three-dimensional buildings results in separated-Gow 
regions fundamentally different from those about two-dimensional buildings. 
Fluid is convected into and out of separated zones for the three-dimensional 
case while turbulent diffusion is the primary mode of fluid exchange for the 
two-dimensional case. 

(2) Mean streamline patterns for flow about buddings have been proposed 
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which provide a basis for further research into wind magnitudes and dispersal 
patterns near buildings. 

(3) Analytical methods have been developed which provide an effective 
prediction of mean velocity, turbulence intensity, mean temperature and 
mean concentration of pollutants in the wake region downwind of buildings. 
Additional research is needed to further refine these methods. 

Acknowledgments 

The authors are indebted to the Nuclear Regulatory Commission, the 
Electric Power Research Institute and the National Aeronautics and Space 
Administration for support which has led to many of the results discussed in 
this paper. 

References 

1 M. J. Colmer, Some full-scale measurements of the flow in the wake of a hanger, ARC- 
CP-1166, 1971. 

2 W. Frost and A.M. Shahabi, A field study of wind over a simulated block building, 
NASA Rep. CR-2804, 1977. 

3 J. Counihan, An experimental investigation of the wake behind a two-dimensional 
block and behind a cube in a simulated boundary layer flow, CERL Lab. Note RD/L/ 
N115/71, 1971. 

4 R. Lemberg, On the wakes behind bluff bodies in a turbulent boundary layer, Uni¬ 
versity of Western Ontario Rep. BLWT-3-73, 1973. 

5 J.A. Peterka and J.E. Cermak, Turbulence in building wakes, 4th Int. Conf. on Wind 
Effects on Buildings and Structures, London, 1975, Cambridge University Press, 
London, 1977, pp. 447—463. 

6 A.C. Hansen and J.E. Cermak, Vortex-containing wakes of surface obstacles, Fluid 
Dynamics and Diffusion Laboratory Rep. CER75-76ACH-JEC16, Colorado State 
University, Fort Collins, CO, 1975. 

7 I.P. Castro and A.G. Robins, The effect of a thick incident boundary layer on the 
flow around a small surface mounted cube, Central Electricity Generating Board 
Rep. R/M/N795, 1975. 

8 H.G.C. Woo, J.A. Peterka and J.E. Cermak, Wind-tunnel measurements in the wakes 
of structures, NASA Contractor Rep. NASA CR-2806, 1977. 

9 J.C.R. Hunt, C.J. Abell, J.A. Peterka and H.G.C. Woo, Kinematical studies of the 
flows around free or surface-mounted obstacles; applying topology to flow visualiza¬ 
tion, J. Fluid Mech., 86 (1978) 179—200. 

10 K.M. Kothari, J.A. Peterka and R.N. Meroney, Stably stratified building wakes, U.S. 
NRC Rep. NUREG/CR-1247, 1979. 

11 R.P. Hosker, Flow and diffusion near obstacles, in D. Randerson (Ed.), Atmospheric 
Science and Power Production, U.S. Department of Energy, Chap. 7, 1984. 

12 M.E. Davies, V.G. Quincey and S.J. Tindall, The near-wake of a tall building block in 
uniform and turbulent flows, in J.E. Cermak (Ed.), Proc. 5th Int, Conf. on Wind 
Engineering, Fort Collins, CO, 1979, Pergamon Press, Oxford, 1980, pp. 289—298. 

13 W.W. Li and R.N. Meroney, Gas dispersion near cubical model building. Part I. Mean 
concentration measurements, J. Wind Eng. Ind. Aerodyn., 12 (1983) 15. 

14 A.C. Hansen, J.A. Peterka and J.E. Cermak, Wind-tunnel measurements in the wake 
of a simple structure in a simulated atmospheric flow, NASA Contractor Rep. NASA- 
CR-2450, 1975. 



38 


15 J.C.R. Hunt and J.P. Smith, A theory of wakes behind buildings and some provisional 
experimental results, Central Electricity Research Laboratory Rep. RD/L/N31/69, 
1969. 

16 J.C.R. Hunt, Further aspects of the theory of wakes behind buildings and comparison 
of theory with experiment, Central Electricity Research Laboratory Rep. RD/L/R/ 
1665,1971. 

17 J.C.R. Hunt, Vortex and momentum wakes behind surface obstacles in turbulent 
boundary layers, private communication. 

18 R.E. Akins, J.A. Peterka and J.E. Cermak, Mean force and moment coefficients for 
buildings in turbulent boundary layers, J. Ind. Aerodyn., 2 (1977) 195—209. 

19 K.M. Kothari, J.A. Peterka and R.N. Meroney, An analytical method for determina¬ 
tion of dispersion in the wake of a power plant, to be published. 

20 K.M. Kothari, J.A. Peterka and R.N. Meroney, The wake and diffusion structure be¬ 
hind a model industrial complex, Nuclear Regulatory Rep. NUREG/CR-1473, 1981. 

21 K.M. Kothari, R.N. Meroney and J.A. Peterka, Nuclear power plant building wake 
effects on atmospheric diffusion: simulation in wind tunnel, Electric Power Research 
Institute Rep. EPRI NP-1891,1981. 



